Engineering of triply entangled states in a single-neutron system 
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We implemented a triply entangled Greenberger-Horne-Zeilinger(GHZ)-like state and coherently 
manipulated the spin, path, and energy degrees of freedom in a single-neutron system. The GHZ-like 
state was analyzed with an inequality derived by Mermin: we determined the four expectation values 
and finally obtained M — 2.558 ± 0.004 ^ 2, which exhibits a clear violation of the noncontextual 
assumption and confirms quantum contextuality. 
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Einstein, Podolsky, and Rosen argued that quantum 
mechanics (QM) is not a complete theory in the sense 
that some results which can be predicted from remote 
measurements are not described by the theory [1]. Bell 
showed that local hidden variable theories (LHVTs) sat- 
isfy some inequalities which are violated by QM [2] , thus 
QM cannot be completed with LHVTs. Experimental 
violations of Bell inequalities on bipartite systems have 
been observed with two- photons (sj, two-ions [3], atom- 
photon systems jsl, and two-hadrons [ij. Moreover, Bell- 
like inequalities can be tested using different degrees of 
freedom of single-particle systems. In this scenario, the 
violation of the inequality does not prove the impossi- 
bility of LHVTs, but the impossibility of noncontextual 
hidden variable theories (NCHVTs) 0,i3. In NCHVTs 
the result of a measurement A is predetermined and is 
not affected by other previous (or simultaneous) measure- 
ments, carried out on the same individual system, of any 
observables mutually commuting with A. Experimen- 
tal violations of Bell- like inequalities with two degrees of 
freedom of single neutrons have been observed [9|. 

Even more apparent conflicts between predictions by 
QM and LHVTs are found by Greenberger, Horne, and 
Zeilinger: entangled states of three or more separated 
systems can lead to predictions nonstatistically in con- 
tradiction to each other TOf|. Indeed, Mermin showed 
that this conflict can be converted into a larger violation 
of a Bell-like inequality between three or more separated 
systems ll| . Experimental test of these inequalities were 
reported, e.g., with the_use of three and four photons 
13 . and four ions [14]: among these, tests of quan- 



tum non-locality on many-particle generalizations of the 
GHZ triplet are particularly appealing [l^. A natural 
question is whether a violation of Mermin-like inequali- 
ties can be observed also on single-particle systems. The 
interest of this violation goes beyond the technical chal- 
lenge of preparing GHZ-like entangled states using three 
degrees of freedom of a single-particle system and the 
capability of measuring the corresponding observables. 
The violation of the Mermin-like inequality is interesting 
in itself since it is more robust to noise than previous 



violations of bipartite Bell-like inequalities and thus em- 
phasizes the conflict between QM and NCHVTs. 

Starting from a demonstration of a violation of Bell- 
like inequality [9|, several neutron optical experiments 
were accomplished using Bell-like states, with entangle- 
ment of two, i.e., the spin and the path, degrees of 
freedom of neutrons [17|. Recently we developed a 
coherent-manipulation method of a neutron's energy, i.e., 
total energy of neutroi is g iven by the sum of kinematic 
and potential energies |l8|. This technique accompanied 
by phase manipulations [19| allows us to add one more 
degree of freedom to be entangled: a triply entangled 
GHZ-like state in a single neutron system is generated 
and manipulated. Here we report the first preparation of 
a GHZ-like state using three degrees of freedom of a sin- 
gle neutron, i.e., two internal degrees of freedom (the spin 
and the energy) and one external one (the path taken by 
the neutron in an interferometer setup), and the first vi- 
olation of a Mermin-like inequality with a single-particle 
system. General descriptions of perfect crystal neutron 
interferometer experiments can be found in a book [2^ . 

In the neutron interferometer experiments accompa- 
nied by two radio-frequency (RF) oscillating fields (see 
FigdJ, the total state consists of the neutron state I^'at) 
and the two RF fields, \a^) and \a^/2) represented by 



coherent states: l^'/of) — 



<j/2} 
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first, the incident neutron was polarized, e.g., to up, de- 
noted by It): all states, corresponding to the spin, path 
and energy degrees of freedom, are represented by the 
north-pole points of the Bloch-sphere in Fig.l. In passing 
through the first plate (the beam splitter) of the inter- 
ferometer, the state describing neutron's path is trans- 
formed into a 50/50 superposition of path-I (|I)) and 
path-II (|II)) states. Therefore, the corresponding state 
lies on the equator of the path Bloch-sphere. In the inter- 
ferometer, a RF spin-flipper operated with frequency to 
was inserted in the path II, where the spin-flip process by 
a time-dependent interaction induces energy transitions 
from the initial energy state |£'o) to states \Eo — huj) by 
photon exchange: up-spin |t) only in path-II was flipped 
to down-spin ||), thus losing energy by hu; [^H. Conse- 



quently, one can generate the state of neutrons in a triply 
entangled GHZ-like state, given by 

1*^^^) = -l={\mmEo)+\i)m)^\Eo'f^)}- (i) 

Note that, in this GHZ-like state, a product state |t 
) i^i |I) (8" \Eo), all states on the north poles, and another 
product state |J,) (g) |II) (X) \Eo — fku), all states on the 
south poles of the Bloch-spheres, are superposed. Here, 
the state of neutrons is characterized by three, i.e., the 
spin, path and energy, degrees of freedom, which are sim- 
ply described by two-level quantum systems such as 



I ^ spin} 
I'i'path) 



{|T>,U>} 
{|i>,|n)} 

= {\Eo-hLj), 



\Eo)}. 



(2) 



In this simple description, all subspaces are effectively 
spanned by orthogonal two-bases. It is noting here that 
the energy subspace (and the momentum subspace for 
the path degree of freedom) is(are), in principle, not a dis- 
crete two-level but has(have) a continuous structure. Ob- 
servable of one subspace commute with that of a different 
subspace, which justifies the derivation of a Mermin-like 
inequality due to NCHVTs. 

The important operations of each degree of freedom in 
the experiments were phase manipulations between each 
of two bases, (a) The spin-phase a was adjusted by a 
magnetic field oriented along the quantisation axis, i.e., 
+z direction, tuned by an 'accelerator' DC coil. In real- 
ity, the change of the Larmor frequency Auj^ results in 
a phase shift a = Aw^Ti, where Ti is the propagation 
time through the 'accelerator' coil, (b) The phase ma- 
nipulation of the path subspace was accomplished with 
the use of an auxiliary phase shifter made of a parallel- 
sided Si plate 5mm in thickness. In this case, the phase 
shift X 'WS'S given by x = NbcXD, with the atom den- 
sity N, the coherent scattering length be, the wavelength 
of the beam A, and the thickness of the plate D. (b) 
There is a suitable method for a phase manipulation of 
the energy degree of freedom, which is known as a zero- 
field precession when two RF flippers (operated at 
a frequency Ur) are set in serial, the former induces the 
energy difference ±hujr until the latter, resulting in the 
phase difference 7 — 2LUrT2, where T2 is the propagation 
time between the flippers. In particular, an experimen- 
tally convenient method to manipulate individually the 
Larmor phase and the zero-field phase 7 was found and 
reported in [l^. 

Since perfect correlations (or anti-correlations) can not 
be observed in real experiments, one should use an in- 
equality in order to clarify peculiarities of the triply en- 
tangled GHZ-like state. Mermin analyzed the GHZ ar- 
gument in detail and derived an inequality suitable for 
experimental tests to distinguish between predictions by 
QM and by LHVTs Assuming a tripartite system 
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FIG. 1: Schematic view of the main experimental setup (not 
to scale) for the preparation and analysis of triply-entangled 
states in a single-neutron system together with Bloch-sphere 
descriptions to depict evolutions of each quantum state, i.e., 
the spin, path and energy degrees of freedom. The experiment 
consist of three stages, (i) Preparation of a triply entangled 
GHZ-like state |*^^^): the state of neutron is the 50/50 su- 
perposition of It) 55 |I) (8 l-Eo) (all states on the north poles) 
and ID ® |II) (g) \Eo — huj) (all states on the south poles of 
the Bloch-spheres). (ii) Manipulation of the relative phases 
followed by projection measurements: the directions of the 
projection measurements P-' are depicted by thick red ar- 
rows in Bloch-spheres. (iii) Detection: numbers of neutrons 
N(x; a; 7) are counted. 



and taking the assumption in the conditionally indepen- 
dent form (represented by the Eq.(5) in [ll|) due to 
NCHVTs instead of LHVTs, one can obtain the border 
for a sum of expectation values of certain product ob- 
servables, which is to be tested in the experiment. The 
sum of expectation values M is defined as 

M = i?[aPa>^]~i?[a^a^a^]~i?[a^<a;]-i?[a^a^a^] (3) 

where . .], cr^, cr|, and (t| represent expectation values, 
and Pauli operators for the two-level systems in the path, 
spin, and energy subspaces, respectively. NCHVTs set 
a strict limit for the maximum possible value of 2. In 
contrast, quantum theory predicts an upper bound of 4: 
any measured value of M that is larger than 2 decides 
in favor of quantum contextuality. A violation of up to 
factor of 2 is expected with a triply entangled GHZ-like 
state 1^*^^^). 

In order to test the Mermin-like inequality, one should 
determine four expectation values for joint measurements 
of three, i.e., spin, path and energy, observables. We de- 
cided to extend the strategy used in the measurement of 
the Peres-Mermin proof of the Kochen-Specker theorem 
[2^ : successive measurements of three degrees of free- 
dom were carried out. Projection operators P^cj)) to the 
state -^{\^i) + e*"^!*/")} (l=path, spin, energy) were 
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FIG. 2: Typical intensity 
modulations obtained by 
varying the path-phase x- 
The phases a and 7, for 
the spin and the energy, 
respectively, are tuned at 0, 
7r/2, TT and 37r/2 in order to 
accomplish projection mea- 
surements of P(0), -P(7r/2), 
P(7r), and P(37r/2) for both: 
(a) 7 = 0, (b) 7 = 7r/2 (c) 
7 = TT (d) 7 = 37r/2 and in 
each 7 setting a was set at 0, 
7r/2, TT and 37r/2 (upper to 
lower panels). 



involved in the experiments: <j\ = P'(0) — P^tt) and 
a'y = P'(7r/2)-P'(37r/2). Note that these projection op- 
erators differ only in phase between the two orthogonal 
states, which was experimentally very convenient since 
one only needs phase manipulations: directions of all pro- 
jection measurements are depicted by thick red arrows 
in Bloch-spheres in Fig.l (they all lie on the equatorial 
plane). In practice, each expectation value was deter- 
mined by a combination of eight count rates in a single 
detector with appropriate phase settings, for instance, 

^ N(0:^:^)-Ni7r:^:^)~N{0:f:^)...+NiO:f:f)-Ni^:^:f) 
Ar(0:f :f )-(-Ar(7r:f :f )+JV(0:if :f )...+Af(0:^:^)+Af(7r:if iSf) 

(4) 

where N{x;a;-f) = {^\PP{x) ' ^"(a) ' P%l)\'^) denotes 
the count rate with the path phase x, the spin phase a, 
and the energy phase 7. 

The experiment was carried out at the perfect-crystal 
neutron-interferometer beam line S18 at the high flux re- 
actor at the Institute Laue Langevin (ILL). A schematic 
view of the main components of the experimental setup 
together with a Bloch-sphere description depicting evo- 
lutions of each degree of freedom is shown in FiglTJ A 
silicon perfect-crystal monochromator was placed in the 
neutron guide to monochromatize the incident neutron 
beam to a mean wave length of Ao=1.92A with the 
monochromaticity AA/Aq « 0.01. The cross section of 
the incident beam was confined to 5 x hmrn?'. Magnetic 
prisms were used to polarize the incident beam vertically, 
before the beam enters a triple-Laue (LLL) interferom- 
eter. The interferometer was adjusted to give the 220 
reflections. A parallel-sided Si plate was used as a phase 
shifter to tune the phase x for the path degree of freedom. 
This phase shifter accompanied by the beam recombina- 
tion in the interferometer enabled to realize a projection 
measurement with the operator PP(x)- 

A fairly uniform magnetic guide field, Bq in +z 
(~20G), was applied around the interferometer by a pair 
of water-cooled Helmholtz coils (not shown in FiglT]). 
The first RF spin-flipper was located in this region, and 



its operational frequency was tuned to 58kHz. The GHZ- 
like state of neutrons l^'^^'^) was generated by turn- 
ing on this RF spin-fiipper. Along the flight path af- 
ter the interferometer, another fairly uniform magnetic 
guide held, B'q (at half strength ~10G), was appfied with 
another pair of water-cooled coils in Helmholtz geome- 
try (also not depicted in FiglT]). The second RF spin- 
flipper, tuned to the operational frequency of 29kHz, 
was placed in this region. This RF spin-fiipper was 
mounted on a common translator together with a DC 
spin-flipper. The translation of the common basis al- 
lows one to tune the phase 7 of the energy degree of 
freedom independently ^1^; for instance, 7 = 0, 7r/2, tt, 
37r/2 resulted in implementing P%0), P%Tr/2), P''{tt), 
and P'^(37r/2). Note that the second RF spin-fiipper, in 
practice, worked as an energy "recombiner" described as 
= ^\E^^^nuj l2){{Eo\ + {Eo-hu)\}. A spin-analyzer 
in the +z direction (a bent magnetically saturated bent 
Co-Ti supermirror) together with a 7r/2 spin-turner en- 
abled the selection of neutrons in xy-plane (normal to 
the quantization axis). An accelerator coil, oriented in 
Bacc+z was used to adjust the spin phase a = 0, 7r/2, tt, 
37r/2 accomplishing projection measurements of P*(0), 
P"(7r/2), P'{tt), and P"(37r/2). 

By tuning the spin phase a, and the energy phase 7 
each at 0, 7r/2, tt and 37r/2, 16 independent path phase 
X scan, i.e. oscillation measurements, was carried out 
for the determination of M . Typical oscillations are de- 
picted in FigO intensities at indicated fines (x = 0, 7r/2, 
TT, 37r/2) were used to determine the related expectation 
values. The contrasts of the oscillations were just below 
70%, which was about the same as those with the empty 
interferometer and showed that all parameters could be 
manipulated effectively. 

Measured intensity oscillations were fitted to sinusoidal 
curves by the least squares method and the four related 
expectation values were extracted. Statistical errors were 
estimated to ±0.001 taking all fit-errors from single mea- 
surement curves into account. One set of measurements, 
consists of 32 oscillation measurements. (We recorded in- 
tensities with and without spin-fiips at each phase shifter 
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TABLE L Experimentally determined expectation values and 
the resulting M value. 



Observables 


X 


Variables 
a 


7 


Values 


p s e 

^ X ^X ^ X 


0,7r 


0,7r 


0,7r 


0.659(2) 


p s e 
^ X ^y^y 


0,7r 


7r/2,37r/2 


7r/2,37r/2 


-0.603(2) 


^ y^ x^ y 


7r/2,37r/2 


0,7r 


7r/2,37r/2 


-0.632(2) 


p s e 
^ y^ X 


7r/2,37r/2 


7r/2,37r/2 


0,7r 


-0.664(2) 




M = 


2.558 ± 


0.004 





Here, as the first reahzation of triple entanglement in a 
single-particle system, the GHZ-like state was generated 
and analyzed. Now we are ready to proceed to investi- 
gate other triply entangled states, for instance, the W- 
state [iJ] which is expected to be generated rather easily 
with a double-loop neutron interferometer 25 1. 



X position, which allowed estimation and correction, if 
necessary, of the path-phase x instability afterwards.) 
We measured 4 sets of such 32-oscillations to reduce sta- 
tistical errors. During the analysis, we noticed that the 
statistical errors here are much smaller than those ob- 
tained in the Bell-like inequality experiments 0. This 
is due to the fact that the points used to determine ex- 
pectation values are in the vicinity of the flat maxima or 
minima, e.g. N{x : a = 0, tt : 7 = 0, tt) around x = 0, tt 
on the solid and dotted lines in Fig [21 which reflects ro- 
bustness of the Mermin-like inequality and led to rather 
small statistical errors. Four measurements were summed 
up as weighted averages and the final value and the er- 
ror were determined. So, the final errors are the sum of 
systematic and statistical errors. (Systematic errors were 
mainly due to the path-phase x instability, i.e., unwanted 
drifts of the x- phase, during the measurement.) We ob- 
tained four expectation values listed in Tab.l together 
with settings of variables and the final M-value. In eval- 
uating the Mermin-like inequality, M was calculated to 
be M = 2.558 ± 0.004. This exhibits a clear violation 
M ^ 2 of the noncontextual border. The reduction from 
the ideal value of 4 is solely due to reduced contrast of 
the interference term from the interferometer, i.e. just 
below 70%. 

Our results with neutrons were obtained with detec- 
tors of more than 99% efficiency. This experiment alone 
will not close all loopholes, e.g., a light-cone loophole is 
still remaining, but such a high efficiency of detectors 
for neutrons will help to consider physics of contextual- 
ity thoroughly. The use of entanglement of the energy 
degree of freedom is not limited to neutrons but easily 
applicable to other quantum systems. Furthermore, a co- 
herent manipulation of energy degree of freedom can be 
extended to create artificial multi-level quantum system, 
e.g., in the order of lO'^, in a single-particle system by ap- 
plying a multiple-frequency energy-manipulation scheme 
in serial. Such a system could be used for quantum in- 
formation processing. 

In summary we have accomplished the demonstration 
of the violation of the Mermin-like inequality with the use 
of three, i.e., the spin, path and energy, degrees of free- 
dom in a single-neutron system. The concept of entangle- 
ment is not limited between spatially-separated systems 
but also generally applicable between degrees of freedom. 
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